„ono m/< ^ 



Mem. S.A.It. Vol. 84, 1 
©SAIt 2008 



Memorie deiia • 




e ^OegliS9* 



The cosmic X-ray background: abundance and 
evolution of hidden black holes 
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Abstract. The growth of supermassive black holes across cosmic time leaves a radiative 
imprint recorded in the X-ray background (XRB). The XRB spectral shape suggests that a 
large population of distant, hidden nuclei must exist, which are now being revealed at higher 
and higher redshifts by the deepest surveys performed by Chandra and XMM. 
Our current understanding of the XRB emission in terms of AGN population synthesis 
models is here reviewed, and the evolutionary path of nuclear accretion and obscuration, as 
emerging from the major X-ray surveys, is investigated. The role of galaxy merging versus 
secular processes in triggering nuclear activity is also discussed in the framework of recent 
galaxy/black hole co-evolutionary scenarios. Finally, the limits of current instrumentation in 
the detection of the most obscured and distant black holes are discussed and some possible 
directions to overcome these limits are presented. 
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1. Introduction 

Fifty years of technological progresses and sci- 
entific investigations since its discovery, have 
revealed the main essence of the cosmic X- 
ray background (XRB), the diffuse X-ray glow 
pervading the sky discovered by Giacconi et 
al. in 1962 in the first X-ray astronomy exper- 
iment. From that first rocket flight to the deep- 
est observations now performed with Chandra 
and XMM, X-ray surveys gained ten orders of 
magnitude in sensitivity. This progress showed 
that, at energies below 10 keV, the XRB glow 
is resolved into hundreds of millions of in- 
dividual X-ray sources distributed across the 
whole Universe: most of them are accreting 
supermassive black holes shining as Active 
Galactic Nuclei (AGN) at different epochs. In 



Send offprint requests to: R.Gilli 



the center of the deepest X-ray field to date, the 
4Ms Chandra Deep Field South (4Ms CDFS; 
Xue et al. 2011), an AGN surface density of 
~10 4 deg~ 2 is observed down to f2-uev~ 5 x 
10~ 17 erg cirT 2 s -1 , whose integrated emission 
can virtually account for most of, if not all, the 
XRB surface brigthness below 10 keV (Xue 
et al., 2011; Moretti et al., 2012). At higher 
energies, however, in the 20-40 keV energy 
range, where the XRB spectrum peaks, only a 
tiny fraction of its emission has been resolved 
into individual sources, mainly because of the 
limitations of high-energy instrumentation. A 
large portion of the XRB emission at 30 keV is 
expected to be produced by the same sources 
observed in the deep X-ray surveys like the 
4Ms CDFS, just by simply extrapolating their 
spectra to higher energies. However, the very 
shape of the XRB spectrum calls for an ad- 
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Fig. 1. Compilation of the measurements of the cosmic X-ray background spectrum in the 0.5 - 400 keV en- 
ergy range. Datapoints with different colors come from different combinations of missions and instruments 
as labeled and referenced (left labels: E<10 keV; right labels: E>10 keV). 



ditional large population of heavily obscured, 
Compton-thick AGN (hereafter CT AGN, N H > 
10 24 crrT 2 ) poorly sampled by most surveys at 
E<10 keV. Despite several efforts, the cosmo- 
logical evolution and luminosity function of 
CT objects remain essentially unknown and 
have to be postulated by AGN population syn- 
thesis models trying to explain the XRB broad 
band spectrum (see next Section). Heavily hid- 
den black holes are then likely the key and 
still missing ingredient to get a complete un- 
derstanding of the cosmic XRB, keeping it one 
of the most fascinating and challenging topics 
in high-energy astrophysics. 

2. The synthesis of the X-ray 

background: Compton-thick AGN 
and the "missing" fraction 

The broad band, 0.5-400 keV XRB spectrum 
has a characteristic shape peaking at E~30 keV. 
Below 10 keV it can be approximated by a 
power law with T~1.4, i.e. it is harder than the 
average spectrum of bright, unobscured QSOs. 
A compilation of XRB measurements includ- 



ing the most recent results by Chandra and 
XMM at E< 10 keV, and BeppoSAX, Swift and 
INTEGRAL at E>10, keV is shown in Fig.l. 
The various measurements generally agree on 
the spectral shape, but possess significant scat- 
ter, at the 20-30% level, in their absolute nor- 
malization (in general, most of the recent XRB 
measures were found to be higher than the clas- 
sic value measured by HEAO-1 in the 80s). 
The origin of this scatter is still debated. Part 
of it, but limited to the measures performed 
in pencil beam surveys, could arise from cos- 
mic variance. Another source of uncertainty 
could be the stray-light which is entering an 
X-ray telescope field of view if not properly 
modeled (Moretti et al., 2012). Finally, >10% 
calibration uncertainties have been observed 
among different instruments on board different, 
or even the same, missions (Tsujimoto et al., 
201 1), which could explain another part of the 
scatter. It is then fair to conclude that the XRB 
absolute flux is known with a ~20% systematic 
uncertainty. 

In 1989 Setti & Woltjer first proposed that 
the hard XRB spectrum could be explained 
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Fig. 2. Fraction of the XRB emission at 20 keV pro- 
duced by CT AGN at different redshifts in the Gilli 
et al. (2007) model. The emission from the whole 
CT population (black line) is splitted into four bins 
of intrinsic 2-10 keV luminosity L"" (color curves): 
most of the "missing" XRB is expected to be pro- 
duced by objects with L"" < 10 44 erg s _1 and z < 1. 



by the superposition of AGN spectra with dif- 
ferent absorption degrees. This prompted the 
flourishing of a series of AGN population 
synthesis models with an increasing degree 
of complexity, including: detailed spectra and 
absorption distributions (Madau et al., 1994; 
Comastri et al., 1995); a careful treatment 
of Compton scattering in heavily obscured 
sources (Wilman & Fabian, 1999; Pompilio 
et al., 2000); the inclusion of iron features 
(Gilli et al., 1999); evolution of the obscured 
AGN fraction with luminosity (Ueda et al., 
2003) and redshift (La Franca et al., 2005; 
Ballantyne et al., 2006; Treister et al., 2009); 
dispersion in the AGN primary spectral slopes 
(Gilli et al., 2007); exploration of the full pa- 
rameter space (e.g. reflection fraction vs CT 
abundance; Treister et al. 2009; Akylas et al. 
2012). These works suggested that 10-30% of 
the XRB emission at 30 keV, depending on the 
model assumptions and on the adopted XRB 
normalization, is produced by CT AGN. Gilli 
et al. (2007) splitted the population of CT AGN 
into a pair of evenly populated categories of 
transmission-dominated (NH~10 24 crrr 2 ) and 
reflection dominated (NH>10 25 crrT 2 ) objects, 




Fig. 3. The distribution of the observed 2-10 keV 
to [Ne V] luminosity ratio (X/NeV) for the [Ne V]- 
selected Type-2 AGN in COSMOS. Filled and 
empty histograms refer to X-ray detections and up- 
per limits, respectively. The filled dot (X/NeV~9, 
arbitrary position on the v-axis) is derived by stack- 
ing X-ray undetected sources. The vertical dashed 
line marks the threshold to select candidate CT AGN 
(X/NeV<15). From Mignoli et al. (2013). 



and built a self-consistent synthesis model in 
which CT AGN are as abundant as less ab- 
sorbed objects. However, this result is far from 
being conclusive: as Comastri et al. (2008) 
showed, varying the number of reflection- 
dominated objects by a factor of four was still 
consistent with the XRB spectrum. This sug- 
gests that the XRB spectrum has limited con- 
straining power and that an accurate estimate 
of the abundance of CT AGN can only be de- 
termined by direct object selection. 

In Fig. 2 it is shown the fraction of the 
XRB emission at 20 keV produced by CT 
AGN (the so-called"missing" XRB) as a func- 
tion of their redshift in the Gilli et al. (2007) 
model. Most of the "missing" XRB is expected 
to be produced by Seyfert-like CT objects at 
z~0.7 (i.e. with intrinsic 2-10 keV luminosity 
L'"'<10 44 erg s~'). A direct sampling of this 
very population is therefore highly desirable, 
and it has been recently attempted by Mignoli 
et al. (2013) who used the [Ne V] 3426 A emis- 
sion line to pinpoint narrow line (type-2) AGN 
among the -7400 galaxies at0.65<z<1.2inthe 
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zCOSMOS spectroscopic survey (Lilly et al., 
2007). Seventytwo type-2 AGN with Seyfert- 
like luminosities and average redshift 0.9 were 
found in the survey area covered by Chandra 
(~0.9deg 2 , the C-COSMOS survey, Elvis et al. 
2009): among them, 46 objects with X-ray to 
NeV ratio lower than 100 were accounted as 
heavily obscured AGN, including 9 CT candi- 
dates with X/NeV< 15. The distribution of the 
X/NeV ratios for this sample is shown in Fig. 3: 
by stacking the X-ray undetected sources we 
derived an average X/NeV ratio of 9, which 
points towards a CT fraction in the whole type- 
2 sample of * 41% (see the detailed simu- 
lations in Mignoli et al. 2013), in line with 
the 50% fraction assumed in the Gilli et al. 
(2007) model. This measurement will be fur- 
ther refined by the forthcoming 2.8Ms Chandra 
COSMOS-Legacy data (PI F. Civano). 

3. Evolution of nuclear obscuration 
and triggering mechanisms 

While an overall self-consistent picture is 
emerging on the CT AGN population from the 
previous estimates, it is recalled that any esti- 
mate of nuclear absorption based on the com- 
parison between some proxy of the intrinsic 
power (e.g. high-ionization lines such as [O III] 
or [Ne V], or mid-infrared luminosity) and the 
observed X-ray emission is highly uncertain, 
and in turn affects the estimate of the space 
density of CT AGN (see e.g. Fiore et al. 2009; 
Vignali et al. 2010; Georgantopoulos et al. 
2011; Alexander et al. 2011). The only un- 
ambiguous way to select bona-fide CT AGN 
is through X-ray spectroscopy, and deep X- 
ray surveys are key to this. As an example, 
several CT AGN, up to z~5, have been re- 
ported in the CDFS based either on Chandra or 
XMM data (Norman et al., 2002; Tozzi et al., 
2006; Georgantopoulos et al., 2007; Comastri 
et al., 2011; Feruglio et al., 2011; Gilli et al., 
2011). These numbers are too small to com- 
pute their space density and evolution, and this 
situation is unlike to change unless new X-ray 
facilities become available (see next Section). 
Nonetheless, X-ray spectroscopy in the deep 
Chandra and XMM surveys is shedding light 
on another very controversial issue, i.e. the 



evolution history of Compton-thin obscura- 
tion. Several works in the past years claimed 
that the fraction of obscured nuclei increases 
with redshift (La Franca et al., 2005; Hasinger, 
2008; Treister et al., 2009). This relation was 
not found by others (Ueda et al., 2003; Gilli 
et al., 2007), leading Gilli et al. (2010) to pro- 
pose that strong selection effects were pro- 
ducing most of it. X-ray spectroscopy in the 
Chandra/XMM deep fields is shedding new 
light on this issue. Iwasawa et al. (2012) an- 
alyzed the XMM spectra of 46 AGN at z > 1.7 
in the XMM-CDFS and found that the fraction 
of luminous (L™f w > 10 44 erg s _1 ) obscured 
QSOs was significantly higher than that ob- 
served in the local Universe by Burlon et al. 
(201 1), and that selection biases could not ex- 
plain this difference. Based on the Chandra 
spectra of a sample of 34 AGN at z > 3 in 
the 4Ms CDFS, Vito et al. (2013) found that 
the increase of the obscured AGN fraction with 
redshift was less significant or even absent at 
lower luminosities. This behaviour has an in- 
teresting qualitative interpretation in a scenario 
involving two modes of accretion/obscuration, 
which could mirror the two modes of star 
formation (bursting vs secular/main sequence) 
observed in the galaxy population (Elbaz et al., 
2011; Rodighiero et al., 2011). Indeed, bright 
QSOs are preferentially found in starburst 
systems, while lower-luminosity objects are 
preferentially seen in main sequence galaxies 
(Rovilos et al., 2012). Furthermore, the frac- 
tion of AGN and star forming systems in merg- 
ers increases with their bolometric luminos- 
ity and "starburstiness", respectively (Treister 
et al., 2012; Kartaltepe et al., 2012). Based on 
all this evidence, Iwasawa et al. (2012) sug- 
gested that the luminosity-dependent increase 
of the obscured AGN fraction with redshift 
could be related to the higher merger rate and 
gas fraction f gas in galaxies at high redshift. 
As often postulated (e.g. Hopkins et al. 2008), 
mergers/interactions could trigger both pow- 
erful nuclear accretion and a starburst with 
compact and chaotic geometry, installing an 
evolutionary sequence in which, after an ini- 
tial obscured phase, the black holes quickly 
(0.01-0.1 Gyr) clears its close environment 
and shines as an unobscured QSO. The higher 
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Fig. 4. Simulated WFXT spectrum of CDFS-202, a 
CT AGN at z=3.7 (Norman et al., 2002; Comastri 
et al., 2011). About 500 net counts are collected in 
400 ksec (the exposure of the WFXT deep survey). 
The iron line is clearly detected and allows an accu- 
rate redshift determination. 



merger rate alone, however, would simply in- 
crease the abundance of the total QSO pop- 
ulation at high-z. It is the higher gas fraction 
that would produce a longer obscured phase at 
high-z and hence a higher fraction of obscured 
QSOs. If lower-luminosity objects are instead 
preferentially hosted by main sequence, non- 
interacting galaxies, the overall system geom- 
etry could be regular. Hence, a static geometry 
in the nuclear regions, as postulated in the clas- 
sic Unified Models, could be valid and the ob- 
served obscuration would only depend on the 
line of sight. In this case, although a larger gas 
reservoir (higher f gas ) is available at high-z, the 
lack of interactions would keep it away from 
the BH: the available gas would then just be 
consumed to form stars. The system geometry 
will then be preserved in time and, in turn, the 
observed fraction of obscured AGN. This ap- 
pealing scenario admittedly relies on a qualita- 
tive interpretation: further observations and de- 
tailed modeling are needed to test it in the next 
few years. 

4. Prospects for the very near and 
mid-term future 

Besides further analysis of Chandra and XMM 
data from completed (CDFS, COSMOS, 
AEGIS, XBootes) and on-going Ms sur- 



veys (COSMOS-Legacy, XXL), a fundamen- 
tal contribution to our understanding of hidden 
black holes is shortly expected by NuSTAR 
(Harrison et al., 2010), on orbit since June 
2012. NuSTAR is carrying the first focus- 
ing optics at E>10 keV and will overtake the 
sensitivity of current high-energy instruments 
by more than two orders of magnitude. The 
NuSTAR surveys in the CDFS and COSMOS 
fields will likely reach sensitivities of a few 
xl0~ 14 erg cur 2 s -1 in the 10-40 keV range. 
These appear obviously very shallow if com- 
pared with the Chandra or XMM data in the 
same fields. It is then unlikely that NuSTAR 
uncovers new X-ray sources in the CDFS and 
COSMOS, but high-energy data will signifi- 
cantly improve the Nh estimate for those al- 
ready known. As far as the demography of CT 
AGN is concerned, the whole NuSTAR pro- 
gram is expected to detect Seyfert-like objects 
(L x <10 44 erg s _1 ) up to z<0.5 (in small num- 
bers though), opening a new territory to test the 
assumptions of XRB synthesis models. 

A real breakthrough in the study of the 
demographics of CT AGN would be possi- 
ble in the next years if concept-missions like 
the Wide Field X-ray Telescope (WFXT; see 
Murray et al. this volume) will come into real- 
ity. WFXT would produce wide-and-deep sur- 
veys in the 0.5-7 keV band with an unprece- 
dented grasp and average resolution across the 
field of view comparable to that of Chandra. 
Unlike standard Wolter-type mirrors, the poly- 
nomial mirrors desing of WFXT is optimized 
to provide a constant PSF with off-axis angle, 
making this instrument a perfect survey ma- 
chine. To give a sense of its capabilities, with 
its current survey strategy WFXT is expected 
to produce the equivalent of ~1000 CDFS sur- 
veys (at the lMs CDFS depth) plus 3000 C- 
COSMOS surveys. Based on the predictions of 
current XRB synthesis models, WFXT is ex- 
pected to return ~ten millions of AGN up to the 
highest redshifts (z > 6): about 5 X 10 5 AGN 
will be detected with more than 400 net counts, 
including 10 4 objects with Nh > 10 23 cm -2 , 
for which accurate X-ray spectra will be ob- 
tained. This spectral sample will allow an un- 
ambiguous determination of how unobscured 
and moderately obscured AGN evolve as a 
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Fig. 5. Minimum bolometric luminosity of unob- 
scured AGN at z > 6 detectable over a given sur- 
vey area. Different colors refer to surveys in differ- 
ent wavebands (see upper-left corner). Shaded areas 
with progressively darker colors show the regions 
where a given survey should detect more than 100, 
1000, and 10000 objects (assuming an exponential 
decline of the space density of AGN at z >3; Brusa 
et al. 2009; Civano et al. 2011. WFXT matches in 
area and sensitivity future optical (HyperSuprime 
and LSST) and near-IR (Euclid) surveys. 



function of redshift, luminosity, column den- 
sity and environment, providing key informa- 
tion on BH/galaxy coevolution theories. In this 
sample, about 300 objects are expected to be 
bona-fi.de CT AGN at z > 1. As shown by 
the simulation in Fig.4, their WFXT spectra 
will reveal directly their CT nature, and mea- 
sure their redshift through the detection of the 
iron Kq- feature. The cosmological evolution of 
the elusive CT AGN population will then be 
probed for the first time. Finally, WFXT is ex- 
pected to detect > 1000 obscured and unob- 
scured AGN at z > 6 (current model predic- 
tions vary by orders of magnitude), hence de- 
termining how the first black holes evolve: as 
shown in Fig. 5, these data will perfectly match 
the other wide-and-deep multiwavelength sur- 
vey data available in a 5-10 years time frame. 
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